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We report on a study of dimuon events produced in pp collisions at ^/s = 1.96 TeV, using 
0.9 fb~^ of data recorded by the DO experiment during 2008. Using information from the inner- 
layer silicon tracking detector, we observe 712±462±942 events in which one or both muons are 
produced in the range 1.6 < r < 10 cm, which is expressed as a fraction (0.40 ± 0.26 ± 0.53)% 
of the total dimuon sample. We therefore see no significant excess of muons produced a few 
centimeters away from the interaction point. 



1 Introduction 

A recent study by the CDF collaboration claims to observe a large sample of dimuon events in 
which one or both muons appear to be produced at large radial distances (> 1.5 cm) from the 
primary interaction point, with observations not consistent with current heavy quark production 
models. We describe a corresponding search at DO for dimuon events in which one or both muons 
are produced at radial distances exceeding 1.6 cm, relative to the primary pp interaction. 



2 The DO Detector 

DO is a general purpose colliding beam detector described in detail elsewhere'^. The most 
important detector components for this analysis are the tracking and muon systems in the 
central region corresponding to pseudorapidity |r/| < 1.0. 

The central tracking detector at DO is contained within a 2 T solenoidal field, and comprises 
the silicon microstrip tracker (SMT) and the central fiber tracker (CFT). The silicon tracker 
consists of six barrel/disk modules, two outer disks, and an inner layer (LO). The LO detector's 
surrounds the beryllium beampipe, and is formed from single-sided sensors staggered at radii of 
1.60 and 1.76 cm to provide 98% 4> coverage, in the range \z\ < 38 cm. 

The muon system in the |r/| < 1 region consists of a combination of proportional drift tubes 
(PDTs) and scintillation counters'^. These are arranged on both sides of an iron toroid of 
thickness 1.09 m, which provides a magnetic field of 1.8 T to aid in muon identification and 
reconstruction. The total thickness traversed by muons exiting the iron is 12.8-14.5 A giving a 
minimum momentum of around 3 GeV/c for that muon topology. 



3 Event Selection 



This analysis uses data collected by the DO experiment between August and December 2008, 
corresponding to a total integrated luminosity of around 0.9 fb~^. No particular trigger require- 
ments are enforced (i.e. events are accepted from all triggers, inclusively). The event selection 
scheme is designed to approximately match the corresponding requirements used by CDF in 
their analysis' 

The dimuon sample is produced by selecting the two highest-pT muons in each event, pro- 
vided that they satisfy the following requirements. Both muons must fulfill standard DO quality 
criteria, be associated with muon-system hits on both sides of the toroid magnet, and satisfy 
Pt > 3 GeV/c and \r]\ < 1.0. Both muons must be matched to central tracks, which originate 
within 1.5 cm of each other along the beam axis. The combined invariant mass of the dimuon 
pair, M(/i/i), must lie in the range 5 < M{pL^) < 80 GeV/c^. To remove contamination from 
cosmic ray muons, opposite-sign dimuon candidates are excluded if their azimuthal separation 
A(f> exceeds 3.135 rad. Both muons must also be detected within ±10 ns of the expected arrival 
time (for beam-produced muons) by muon scintillators on either side of the toroid. Finally, 
muons are excluded if they do not pass through the active geometrical limits of the LO detector. 

Two subsets of the dimuon sample are defined using information from the SMT. Muons are 
tagged as "loose" if they have three or more hits in the silicon tracker, and tagged as "tight" 
if they are also associated with a LO hit. A loose (tight) event is one in which both muons are 
loose (tight). Events which are loose but not tight therefore contain at least one muon without 
a LO hit, either due to hit inefficiencies, or production beyond LO. The loose selection has a 
radial range up to 10 cm, determined by examining the hit coordinates of muons accepted in 
this sample. 

The number of events containing one or both muons produced beyond LO, A^(excess), can be 
expressed in terms of the number of loose and tight events, and the dimuon selection efficiency 
of the tight sample relative to the loose, St/l^ 

iV(excess) = iV(loose) - i^i^^li^l . (1) 



4 Tight/Loose Muon Detection Efficiencies 



To measure the relative efficiency of the tight and loose selection requirements, a test sample of 
J/tp ^ fi~^fi~ candidates is selected. The same 0.9 fb^^ data sample described above is used to 
reconstruct J/V' candidates from opposite-charge muons originating at a common vertex, and 
with invariant mass 2.95 < M{n^) < 3.2 GeV/c^. All other requirements applied to the signal 
sample are also applied to the J/tjj muons. The test sample is limited to muons produced within 
LO by enforcing the transverse decay length criterion \Lxy{J < 1.6 cm, which is the case for 
> 99.89% of the selected candidates. 

The resulting test sample contains 119 276 J/^l^ candidates, of which (91.5 it 0.6)% 

are signal J/V' events, as determined by fitting the dimuon invariant mass distribution with an 
appropriate model. From this sample of 238 552 muons, 228 569 satisfy the loose SMT criteria, 
and 210 026 also satisfy the tight (LO) requirement. The mean single-muon efficiency, £t/l{i^)i is 
therefore calculated to be 0.9189 it 0.0006; the equivalent dimuon efficiency is the square of this 
value, Erpji^ = 0.8443 ± 0.0008. Here the uncertainties are purely statistical, and are determined 
using the standard binomial expression. 

The efficiency £t/l{iA can depend on the kinematic and geometrical properties of the muon. 
Figure 1(a) shows the relative tight /loose selection efficiency as a function of transverse momen- 
tum, determined by calculating (tight )/A^ (loose) separately in each px bin. The distribution 




(a) pT-dependent tight/loose selection efficiency, (b) Two-dimensional (2, (^)-dependent tight/loose 
showing the result of a fit to a linear function selection efficiency. The low bins in z are associ- 
(solid line), and the mean value (dashed line). ated with the ~1 mm gaps between the sensors. 



Figure 1: The relative tight/loose selection efficiency for muons from J/i/; decays, as a function of their transverse 

momentum pt and two-dimensional geometry (z, 0). 



suggests that higher-py muons may have higher LO efficiencies, at the ~1% level, which is 
quantified by fitting the resulting distribution to a first degree polynomial, yielding the param- 
eterization e(pr) = a 6 • (pT - 9), where a = 0.9214 ± 0.0011, h = 0.0007 ± 0.0003, and the 
transverse momentum is in units of GeV/c. The fit converges with a minimum \^ of 31, for 22 
degrees of freedom. 

A similar calculation is carried out for the muon pseudorapidity 77. In this case, small (±3%) 
and uncorrelated variations in efficiency are observed in different bins, and are taken into account 
on a bin-by-bin basis as described later. Finally, the efficiency is determined as a function of 



the muon {z,(j)) coordinate, as it passes through the LO detector, as shown in Fig. 1(b), This 
accounts for the effects of varying sensor performance. 

The relative efficiency of loose and tight requirements for dimuon events is then taken as 
the product of the efficiencies for the individual muons, where the individual muon efficiency is 
considered as the product of kinematic and geometrical contributions: 

ei^Jii) = e(z^„ <A^J • ^(r/^J • ^(p^^) . (2) 



The geometrical efficiency £{z^-,(j)^.) is taken directly from the histogram shown in Fig. 1(b) 
The normalized pseudorapidity dependence J'^{r]^^) is taken directly from the appropriate bin of 
the e(/ii) vs. r/ histogram. The normalized transverse momentum dependence is parameterized 
as: 



a + b-{pfp-9) 



= (3) 

with a and b fixed at their central values, and transverse momentum in units of GeV/c. For 
muons with transverse momenta larger than 15 GeV/c, the pT correction is fixed to unity. 

5 Results 

In total, the signal sample contains 204177 dimuon events, of which 177535 satisfy the loose 
silicon hit requirements. Of these, 149 161 also pass the tight SMT requirement that both muons 
have a hit in LO. The determination of A^(excess) proceeds on an event-by-event basis: 

Ar(tight) 

A^(excess) = iV"*^^ (loose) - - — (4) 

i=l ^T/L 



Here iV°'^^ (loose) is the observed number of events in the loose sample. The summation is over 
all tight events, and the efficiency for an event is the product of the two muon efficiencies given 
by Eq. Using this method, the total number, and relative fraction, of events with one or 
both muons produced outside LO is found to be: 



This value is significantly smaller than the corresponding fraction of 12% reported by CDF. 

The uncertainty on the efficiency is determined using ensemble tests, in which the full event- 
by-event calculation of Eq. ^ is repeated 1000 times, with the constituent efficiency factors 
in Eq. ^ allowed to vary pseudo-randomly with an appropriate mean and standard deviation. 
The resulting uncertainty ±0.0024 is then translated into an uncertainty on A^(excess). 

6 Systematic Uncertainties and Consistency Checks 

A systematic uncertainty is assigned to the efficiency calculation by repeating it with different 
binning schemes for the z, and rj histograms, and also with the pT factor removed from the 
parameterization. These sources combine to give a total systematic uncertainty on the fractional 
excess of ±0.53%. 

The study is repeated using only those events which are selected by a dedicated dimuon 
trigger in which the muons are not required to be matched to central tracks. The results are 
consistent with the inclusive trigger selection, with A^(excess) comprising a fraction 0.32% of 
the loose sample (396 / 123,572). 

To check for possible enhancements of A^(excess) in different kinematic regions, it is deter- 
mined separately in bins of M(////), (/)(/x), r?(/i) and pj-ifJ-), and no significant departure from 
zero is observed in any particular region. 

7 Conclusions 

Using 0.9 fb~^ of integrated luminosity, a sample of dimuon events is examined with similar 
event characteristics as those selected in the CDF multi-muon analysis We measure the 
fraction of events in which one or both muons are produced in the range 1.6 < r < 10 cm to 
be (0.40 ± 0.26 ± 0.53)%, significantly smaller than the ~12% observed by CDF. These events 
are expected to have contributions from K and vr decays in flight, as well as residual cosmic 
ray contamination, and further studies are in progress to quantify the composition. Additional 
examination of the dimuon properties indicate no significant excess in any particular kinematic 
regions in rj, (p, pT or M{fifj,). 
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A^(excess) 
iV(excess)/iV°^' (loose) 



712 ± 462 (stat.) ± 942 (syst.) 
(0.40 ±0.26 ±0.53) %. 
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